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' Abstract 

> 

^SJ ' In this paper we study the dynamics of the constrained n-dimensional 

rigid body (the Suslov problem). We give a review of known integrable 
■ cases in three dimensions and present their higher dimensional general- 

\l ' izations. 

o 

cn 

^ ' 1 The Suslov problem 

C^H' The equations of nonholonomic systems are not Hamiltonian. They are Hamil- 

tonian with respect to the ahuost-Poisson brackets [5]. This is the reason that 
the integration theory of constrained mechanical systems is less developed then 
for unconstrained ones. However, in some solvable nonholonomic systems with 
an invariant measure, the phase space is foliated by invariant tori, placing these 
systems together with integrable Hamiltonian systems (see [HI ITI 1^ 1141 IH]). 
In this paper we are interested in the integrability of multidimensional gen- 
^ , eralizations of the Suslov nonholonomic rigid body problem. 

' Consider an n-dimensional rigid body motion around the fixed point O = 

(0, 0, ... , 0) in R". The configuration space is the Lie group SO{n). The matrix 
g S SO(n) maps the orthonormal frame . . . , En fixed in the space to the 
frame ei — g ■ Ei, . . . ,en — g ■ En fixed in the body, ei = (1, 0, . . . , 0)', . . . , e„ = 
(0, ... ,0, 1)*. For a motion g(t) £ SO(n), the angular velocity and momentum 
(in body coordinates) are fl{t) — g^^ ■ g{t) and M — J{Vl) respectively. J : 
so{n) —>■ so(n)* is the inertia tensor and has the form: M = J(f2) = IQ + QI, 
where / is symmetric nxn matrix called mass tensor (see [Hj). Here we identified 
so{n) and so{n)* by the Killing scalar product. 

Suppose that additional left-invariant constraints 

{a\n)^0, z = l,...,r (1) 
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and potential force V(r), where F is a unit vector fixed in the space, are given. 
Let Z) = {ri G so(ri), (a% fi) = 0, i ~ l,...,r} and V be the corresponding 
left-invariant distribution of TSO{n). The distribution V is nonintegrable (i.e., 
constraints are nonholonomic) if and only if D is not a subalgebra of so{n) . 

The equations of motion of a rigid body subjected to the constraints ^ and 
the potential force V^(r) are 

dV ^ 

M = [M,r!] + — Ar + ^A,a\ T = -f! • T, (2) 

i—l 

where Lagrange multipliers Ai are determined from the constraints To 
describe a motion g{t) on SO{n) one should solve an additional equation 17 = 
■ 9- 

For 71 = 3, r = 1, the equations (^3), @ are equations of classical Suslov 
problem (see j^l^^D)- Therefore, the above system can be considered as a 
natural multidimensional generalization of the Suslov problem. 

If D is an eigenspace of J then the equations Q take the form 

dV 

J(r!)=pr^[J(r!),12]+pr^(— AF), T = ^flT, QeD (3) 

Contrary to the general case system (O preserves the standard measure 
dndr onM ^ D{n} X S'"-Hr}, 5"-^ = {r e K", (F,F) = l}. The energy 
E = i(Jri, il) + V{T) is always an integral of Q, ®- If, besides the Hamilto- 
nian function, we have dimD + n — 4 additional independent integrals, by the 
Jacobi theorem the solutions of Q can be found by quadratures. Under the 
compactness conditions, the phase space A4 is foliated by invariant tori with 
quasi-periodic dynamic (after appropriate change of time), similarly as in the 
Liouville theorem |14l . 

Therefore, it is natural to call the system completely integrable if it can be 
integrated by the Jacobi theorem; or more generally (see HOI)- if the trajectories 
of the system belong to invariant tori with the dynamic of the form 

^^='i>(^l,...,^.)"-"^'=-<i>(^l,...,^.)- 

The flow of in the new time dr — ^^^dt is quasi-periodic: dipi/dr = uji. 
Moreover, for almost all frequencies uji, . . . ,ujk, by smooth change of variables 
di = Si(ip) equations Q can be reduced to the form di — Vii — ti-ii/n, i = 
1, . . . , fc, where H denotes the total measure of T'' (see pj). 

Note that the integrability of equations ^ without the constraints (for ex- 
ample the multidimensional Euler T^, Lagrange 1^ and Klebsh cases) im- 
plies non- commutative integrability of unconstrained system. The phase space 
T*SO{n) is foliated on d < dim so(n) dimensional invariant isotropic tori with 
quasi-periodic dynamics (see |17[ fll 15]). In the nonholonomic case there is no 
Poisson structure. In order to precisely describe dynamic on the whole phase 
space V we have to solve kinematic equation g^^{t) ■ g{t) = fl{t). This problem, 
for the Fedorov-Kozlov integrable case, is studied by Zenkov and Bloch |21j . 
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In the next section we give a review of the known integrable cases in three 
dimension. Then we present their natural higher dimensional generalizations. 
It appears that these systems are also completely integrable (section 3). 

2 Three— dimensional integrable cases 

For n = 3 we can identify the Lie algebra (so(3), [•,•]) with (M"^, x) via fl — 
(r^i, r22, i^a) = (— i^23, i^i3, — ^^12)- Then the equations of the motion (O be- 
come; 

dV 

.m ^ .m xn + r X — + Xa, T = Txn, (a,r2) = o, (5) 

Without loss of generality we can take a = (0,0,1), i.e., the constraint is 
fls = — fii2 = 0. Then D = M.'^{Ui,^2}- Geometrically, this means that only 
infinitesimal rotations in the planes ei A 63 and 62 A 63 are allowed. 

If F = 0, then the equations © form a closed system in variables fli, fl2- If 
a is not an eigen vector of J then there is an asymptotic line / in D. The angular 
velocity is moving along the energy ellipse £ — {E{n) = h} asymptotically to 
intersections points W- and w+ of £ with / (see Suslov QHI)- That system, 
obviously, has not an invariant measure. When a is an eigenvector of J, then 
system (jSJ preserves the measure dfidT and ior V = the solutions are very 
simple fl — const. 

Thus, further we assume that a is an eigenvector of the inertia tensor J. 
We can take J = diag(Ji, J2, J3). For the integrability, we need one additional 
integral, independent of the Hamiltonian function. 

The known integrable cases are: 

• The Lagrange case (see Kozlov QJ), when Ji = J2 and V — The 
integral is F = {Jn,T). 

• The Kharlamova case [3 U] where V{T) = BiTi + B2T2. Then F = 

JlillBi + J2^2B2- 

• The Klebsh-Tisserand case (Kozlov [H U) with ^(r) = f(Jr,r) and 
F = ^{jn,,m) - i(^r,r}, where A = eJ-idet J. 

Let us note that in the first case the Lagrange multiplier is A = 0, and this 
is actually an invariant subsystem of the usual Lagrange rigid body system. 
Also, Ziglin proved that if ^(r) = i33r3 than system Q has a complementary 
meromorphic first integral in the complexified phase space Mc = C^jfii, ^2} x 
5c{r} only in the Lagrange case Ji = J2 [221 ■ 

The Kharlamova and Klebsh-Tisserand cases can be uniquely described and 
generalized in the following way (see |7]). If Ji 7^ J2 then the equations lO 
of the Suslov problem are integrable for the potentials: V{T) — i;i(ri,r2 -f 
r§) -|- V2{T2, -|- r§) where vi and V2 are arbitrary functions of two variables. 
The corresponding third integral is: F = i(jrj, ,/f2) + J2Ui(ri,r^ + Tj) + 
JiV2{T2,Tl+Tl). 
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In addition, there are integrable cases of the Suslov problem, which include 
the gyroscopic force with the momentum eF x fi added to the right hand side 
of (O. Since the gyroscopic force is conservative, the Hamiltonian remains to 
be the first integral. As well as the usual Suslov problem, the system with the 
gyroscopic force is integrable for the potentials V = 0, V{r) — BiTi + B2T2 
and V{T) = \{ai-tl + 027! + a^ll) (see 13). 

3 Multidimensional integrable cases 

Without potential force, the equations (0, lO represent a closed system in vari- 
able {Euler-Poincare-Suslov equations). We can consider Euler-Poincare- 
Suslov equations on other Lie algebras as well. They are reduced equations 
of nonholonomic systems on Lie groups with left-invariant kinetic energy and 
left-invariant constraints (see |14[ \V6\). Integrable examples can be found in 

miiniiii]. 

Since in the three-dimensional case, only infinitesimal rotations in the planes 
ei A 63 and 62 A 63 are allowed, Fedorov and Kozlov suggested that multidimen- 
sional analogue of Suslov's conditions can be define in the following way: only 
infinitesimal rotations in the planes ei A e„, . . . ,e„_i A e„, i.e., in the planes 
containing the vector e„ are allowed Therefore, we have the following con- 
straints imposed to the components of the angular velocity matrix: 

^^J=0, i,i<n-l. (6) 

Then 

M = M"-H^^ln, ^2n. . . . , ^n-l,^} X S^'^ {T} . (7) 

Really, for y = 0, the system has the same behavior as three-dimensional 
system. If D is an eigen space of the inertia operator J then the solutions are 
simply f2 = const. In general, all trajectories lying on the energy ellipsoid 8 = 
i(Jf2,J7) = h are double asymptotic: they tend to the diametrically opposed 
points w+ and W- as t ^ icx) (see H). 

With the choice of constraints © we shall see that three-dimensional inte- 
grable cases with invariant measure have their multidimensional analogues. 

Suppose that the mass tensor is diagonal / = diag(/i, /2, . . . , In)- Then D 
is an invariant subspace of J. The orthogonal complement of D in so{n) is 
so{n — 1). Since (so(n), so{n — 1)) is a symmetric pair {[soin — 1), soin — 1)] C 
so{n — 1), [so{n — 1),!?] C -D, \D,D\ C so(n — 1)) the equations (jSJ take the 
following simple form: 

dV 

or coordinately: 

dV dV 
OTi dTn 
F^ F^ri^.^ , i — 1, 2, . . . , 7i 1, 

Fri = TiVLin + T2^2n + ■ ■ ■ + F.,i_iil„_i^„. (8) 
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3.1 The Lagrange case 

Let the mass tensor be of the form / = diag(/i, . . . , /i, J„), the rigid body be 
placed in a homogeneous gravitational field and the center of mass be on the axes 
of the dynamical symmetry e„. Then the potential is of the form V{T) = i?„r„ 
and system ((S)) can be seen as an invariant subsystem of the unconstrained 
n-dimensional rigid body motion: 

{li + = {It - Ij)^ ^ik^kj =0, i,j <n-l, 

k 

[h + In)^in — (h — In) ^Ik^kn — BrJ^m 

k 

i^i = -^^ik^k, i=l,...,n. (9) 

k 

The Hamiltonian equations ^ represent the natural generalization of the three- 
dimensional Lagrange rigid body system. Their integrability is proved by Bel- 
jaev However, this does not give us directly the integrability on the invariant 
subspace My = 17^ = 0, i, j < n — 1. In order to prove the integrability of 
the Suslov problem we shall use the symplectic reduction, instead of verifying 
independence of Beljaev's integrals on the invariant subspace. 

Consider the Lagrange rigid body system on the whole phase space T*SO{n). 
It is invariant with respect to the ri(7/it Hamiltonian action of the group S0{n—1) 
(rotations which fix symmetry axes e„). The corresponding moment map $ : 
T*SO{n) so{n)* is given by functions Af^j, i, j < n — 1 considered as the left 
invariant functions on T* SO{n). One can easily prove: 

Lemma 3.1 The reduced system on $"i(0)/5O(n - 1) = T* {SO{n) / SO{n ~ 
1)) = r*5"^^ is the spherical pendulum. 

The spherical pendulum is a completely integrable system. Therefore, the 
invariant submanifold $^^(0) C T*SO{n) of the Lagrange system, as well as 
the rest of the phase space, is almost everywhere foliated by invariant tori with 
quasi-periodic dynamics. (More about the relationships between integrability 
and reductions one can find in The Suslov problem on 2? C TSO{n), 

after Legendre transformation, coincides with the Lagrange system on $~^(0). 
Thus we obtain: 

Theorem 3.1 The Suslov nonholonomic rigid body problem with constraints 
0), the mass tensor of the form I = diag(/i, . . . , /i, /„) and the center of the 
mass on the axes of the dynamical symmetry e„, in a presence of a homogeneous 
gravitational field (the multidimensional Lagrange case) is completely integrable. 

Note that the unconstrained system on T*SO(n) is invariant with respect 
to the left Hamiltonian action of the group SO{n — 1) (rotations which fix F), 
for any inertia tensor J and potential V{T). In this sense, the Lagrange system 
admits two different SO{n — 1) reductions. For n — 3 this corresponds to the 
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fact that two Eulerian angles, the angle of precession tp and the angle of pure 
rotations if, are cyclic variables. Elimination of the variable with the zero 
value of the conjugate momentum p^p = gives the spherical pendulum on 5^ 
(for example, see ^). 



3.2 The Kharlamova case 

The natural multidimensional generalization of the Kharlamova case can be 
defined as the constrained rigid body motion © in a gravitational field with 
the center of mass laying in the plane spanned by ei, 62, . . . , e„_i. Then the 
potential force is of the form V = BiTi + B2T2 + . . . + B„_ir„_i and equations 
(|HJ) become: 

r^ri^n ^ ^ 1,2,...,?^ 1, 

Tri = Tl^^ln + ^2^2n + • • • + Tn-l^n-l,n- (10) 

Note that the singular points of H10|l are given with r„ = r„ = 0. Also, 
equations (|10() have integrals 

F,,(f7) = ^^a,, ^^n,^, z, J = 1, 2, . . . ,n - 1. 

Bi Uj 

Therefore, it is natural to introduce new coordinates wi, . . . , 71, . . . , 7„ 

in the following way: 

— -f^ "In, 

1 1 + o ^1 + -^'i O • O 1 

= ihn S2i„, z = 2,...,n-l, 

Ui til 

7i = 5 — r^H El, 1 = 2, 

7n - r„. (11) 

Then system (|10|l get the simple form 

-^ = 0, i = 2,...,n-l, 

37=7n^i, z = 1,2, . . . ,n - 1, 
at 

Tn _ / Bi \ _ m;;^ / B, 

dt ~ 



i=2 
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Now, we can easily integrate the above system. We have L02 = W2 , • ■ • , = 
uJ^_x- By substitution of the first equation wi — 7„ to the equations with 7i-s 
we get the trajectory as a function of the variable ijJ\. 

7,; - 7O = c^fK - c^°), i = 2,...,n-l, (13) 

where 7°, . . . , 7jJ_i are initial conditions. 
Since (F, F) = 1, we have 
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where P4 is a four degree polynomial depending on initial conditions as param- 
eters. Finally, to get the trajectory as a function of the time t one should solve 
the elliptic integral 



t-tn = 



,0 ±y^Pi{uJi) 

The curve (71 (wi), . . . ,7„_i(a;i)) lies within the ellipsoid 

^f+E rlx 



h+In 



i=2 



On E we have 7n = F„ = 0. Therefore the variable uji varies in the interval 
[^11^2] where ^1 and ^2 are adjacent roots of the polynomial P4 between which 
it takes positive values and (71(^4), . . . ,7„_i(^i)) G S. Thus, the trajectory is 
either asymptotic or periodic with period 



T = 2 



We can summarize the above considerations in the following theorem 

Theorem 3.2 The phase space of the multidimensional Kharlamova system 
is almost everywhere foliated by invariant circles. The periodic motions can be 
expressed as elliptic functions of time t. 



3.3 The Klebsh— Tisserand case 

By multidimensional Klebsh-Tisserand case we mean the nonholonomic rigid 
body motion ^ in the presence of the quadratic potential V{T) = \{BiT\ + 
... + i?„F2): 

{h + In)^in = {Bi — i?„)FiF„, 

F^ F^^r2^^, i 1,2,...,?^ 1, 

Tri = rii7i„ + T2^2n + . . . + r„_iil„_i „. (14) 
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The system (|14|) is similar to the Fedorov-Kozlov integrable case [S]. We 
have integrals 

F,{n, r) = (B, - s„)r2 + (/, + i = i, 2, . . . , n - 1. 

Note that the Hamiltonian function can be expressed in terms of FiS: E = 
i(nB„ + Fi + ... + F„_i). 

The phase space (|7|l of multidimensional Klebsh-Tisserand system is foliated 
by invariant varieties 

7; = {Fi = Ci, . . . , Fn-l = c„}. 

Suppose that Bi > Bn, i = 1, . . . ,n — 1. Then the functions Fi are positive 
and S^. = {Fi = a} are circles in the planes R^{r2i„, F^}, i = 1, 2, . . . , n — 1. 
The invariant surface 7^ is a two-covering of 

sl^x...xSl_^n{rl + ... + rl_,<i} 



determined by the function F„ — ±y 1 — (F^ + . . . + r^_^). The branching 

points of the covering correspond to the zeros of F^. In particular, if constants 
Ci satisfy inequality 

ci , , c„_i 



Bi — Bn Bn-i — B, 



< 1, (15) 



then F„ 7^ and 7^ is the union of two disjoint (n — 1) -dimensional tori (or less 
dimensional tori if some of c.—s are equal to zero) . 

The motion on 7^ can be solved in term of the new time dr = F„dt. Define 
the angular variables (pi on circles S^. by putting 



Using (|14|l . we obtain the following quasi-periodic dynamics on %: 



— =LU, = J , z 1,2, . . . ,n - 1. 

dr \l Ii + In 

The frequencies Ui depend only of and Bi. If the trajectories are periodic 
on one torus, they are periodic on the rest of the tori as well. Also, in the 
original time, system (|14l) takes the form ^ (pi ^ uJi/T~-^. 

Theorem 3.3 //ci, . . . , c„_i satisfy the inequality |J5|) then the invariant vari- 
ety Tc is the disjoint union of two invariant tori. The multidimensional Klebsh- 
Tisserand system \14^ on Tc is completely integrable. 
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